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THERMAL CONTRACTION OF MERCURY 
Han-Shou Liu 


ABSTRACT 

The thermal conduction problem for Mercury with 
radioactive material through a horizontal layer is solved 
in terms of error functions; and the resulting secular 
decrease of moment of inertia about the axis of rotation 
due to thermal contraction is estimated. It is shown that 
the initiation of the tidal evolution on Mercury was com- 
pletely balanced with the effect of the thermal contraction 
during the epoch of cooling in a period of about 48 mil- 
lion years. 
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THERMAL CONTRACTION OF MERCURY 


INTRODUCTION 

Liu (Liu, 1968) has suggested that the trapping of the rotation of the planet 
Mercury was originally affected by the thermal expansion or contraction of its 
figure during solidification. According to this hypothesis, Mercury was origi- 
nally at a very high temperature, and different parts have since cooled bv dif- 
ferent amounts. Two bulges in a state of thermal strain were thus set up on the 
surface. This suggestion, as contained in Liu’s paper (Liu, 1969), was based 
on results of quantitative calculations. In the present paper, the derivation of 
the results on the underlying physical problem is given. On the basis of the 
analysis in the course of this study, it is shown that the thermal process on the 
planet Mercury during the epoch of solidification in cooling is quite adequate to 
account for the origin of trapping. 


COOLJNG 

Apart from certain complications the radiation from the surface of Mercury 
when in a liquid state could have disposed of the internal heat in a very short 
time (Liu, 1969). Liquids in general contract and become denser as they cool, 
and the material cooled by radiation at the surface sank through the hotter 
liquid below, thus maintaining irregular convection currents and a continual 
supply of heat to the surface. In this stage, when the surface temperature was 
high, the loss of heat by radiation was so rapid that radioactivity could not be 
supposed to have delayed the formation of the crust. But as soon as the crust 
of Mercury became so stiff as to stop convection radioactivity would become a 
controlling influence. The equation of heat conduction is (Jeffreys, 1959) 
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dt 
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where T is the temperature, p the density, c the specific heat, k the thermal 
conductivity and A the rate of generation of heat per unit volume. If the axis of 
x is taken vertically downwards, the derivatives with regard to the coordinates 
y and z can be neglected. Further, if the heat conductivity is assumed to be 
uniform, the condition becomes 
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and 


-I - h 2 — — = — when 0 < x < i 
3t Bx 2 P c 


-I-h 2 —j = 0 when x > 
dt ox 2 


( 3 ) 


where 



pc 


and i is the depth of the radioactive layer of Mercury. 
Since 


= A(x --t) 2 
2k 


satisfies Equation (2), one obtains 
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T = u + — when x > 

2k 


The value of u, being specified for all values from -<» to +oo f is 


u 



0(x + 2qh v' t ) 
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which will necessarily satisfy the following equation (Fourier, 195C) 


du . 2 d 2 u 


VI | 

dT 


dx 


The boundary conditions are such that at t = 0 


u = fix + 


S - A L ' 2 ~ ( x - l ) j when x i 
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and 


u = fix + S - — — when x > 'l 
2k 


By specifying (x) = <£(-x), one finds 


T = fix + (S - 4—^1 erf + kl 
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where 


. A (l - x) when x < i 
2k 


and q = 0 when x > . 

In Equation (4) the error function is defined by 
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For A = 7.3 = 10 13 cal. cm. -3 sec. -1 , c = 0.2 cal. g. 1 deg. 1 , k = 0.005 cal. cm 
sec." 1 deg." 1 , £ = 2.5 x 10 6 cm., and S = 1400 deg., then 


and 


where 
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h 2 = 0.005 cm. 2 sec. -1 

L = 944 deg. 

Q = 456 deg. 

M = 3 x 10 " 5 deg. cm." 1 

SECULAR DECREASE IN MOMENT OF INERTIA 

Throughout the process of cooling, a certain amount of mechanical adjust- 
ment must have occurred in the shape of Mercury. Consider a shell of internal 
radius r and external radius r + dr. The temperature, which is defined by 
Equation (5), is a function of t and r. The raise of temperature in a definite time 
is of course a very small negative value. The density of the shell will change at 
the same time from the initial value p to p (1 - 3aT) where a is the coefficient 
of the linear thermal expansion. Let the inner radius change to r (1 + £). The 
external radius will become 


r(l + /3) + dr 
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Hence the mass of the shell after the change of temperature is 
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But the mass is unaltered. Thus the equation of continuity is 


2/3 + — (r/3) - 3aT = 0. 
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Given the cooling, this is the equation to determine the value of subject to the 
boundary condition that P is zero at the center of Mercury. Therefore for any 
shell P is determined by the changes of temperature within that shell. If the 
shell simply expanded independently of the interior, the radius would increase 
by a.Tr instead of by r/3, so that the excess 


r(/5 - cxT) = Sr 


is due to stretching. From Equation (7), the stretching is given by 
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The changes must take place in a short time dt, then 



Let x be the depth of a point below the surface, then r = r Q - / where r r is 
the radius of Mercury. As the temperature change only extend downwards through 
a small fraction of the radius of Mercury, (x/r 0 ) 2 and terms of higher order can 
be neglected. Thus 
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Now, the value a can be expressed as 
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e, + € 
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where e ! and are two constants. From Equations (5), (G), (7) and (8), one ob- 
tains the amount of stretching at the surface: 
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For ^ = 7x 10" 6 deg." 1 , e 2 ~ 2.4 x 10 deg. -2 , and t= 1.5 x 10 'sec., 
Equation (9) gives 


8 = - 2. 8 x 10" 3 . 


From this the increase of the radius of Mercury by thermal expansion dur- 
ing the period of cooling of the shell is at once found to be about -7 km. The 
rate of increase of the radius is, then, about-4.7 x 10 _1 ° cm. sec. -1 . Therefore 
the average rate of change of the moment of inertia about the axis of rotation 
during the period of cooling is approximately 
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where M m is the mass of Mercury. 
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AXIAL ROTATION 


During the epoch of solidification in cooling, Mercury would rotate about its 
polar axis with a time-dependent inertial tensor. The rotation of Mercury was 


_d_ 

dt 


d(f + cpy 


'(t) 


dt 


= - N 


(H) 


where N is the couple due to solar tides, f the true anomaly and <P the angle of 
orientation of Mercury relative to the sun. Equation (11) can be rewritten in the 
following form, 


where 
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and a is the major semiaxis, e the orbital eccentricity, G the gravitational 
constant, k 2 the Love number, e the phase lag of the tides and M s the mass of 
the Sun. 

For e = 0.2, G = 6.7 x 10“ H dyn. cm 2 g; 2 , k 2 = 0.02, 2e = 0.005, p = 5.0 
g. cm.- 3 , M /M = 6.0 x io 6 , r /a = 4.1 x io ~ 5 and r 0 = 2.42 x io 8 cm., the 

s m 0 

average value of N/C t ) over a period of revolution is 



= 0 (IO* 23 ) 


rad. sec, -2 . 
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The calculation also includes the reversal of the phase lag during perihelion 
passage. 

At the 3:2 resonance condition, d(f + p )/ dt = (3/2) • n = 1.8 x 10 "" rad. 
sec' 1 . where n is the orbital mean motion. From Equation (10), one obtains 



-0(10 23 ) rad. sec.” 2 . 


Therefore, it is seen that the initiation of the tidal evolution was completely 
balanced with the effect of the thermal contraction during the epoch of cooling in 
a period of about 48 million years. 


CONCLUDING REMARKS 

In the review of the hypothesis of Darwin for the Earth, it should be noted 
that in the initial stage the thermal contraction of Mercury may be more effec- 
tive rotationally than the solar tides because the loss of heat at that stage was 
quite exceptionally great. Therefore the estimation of If/C(t) may be regarded 
as a lower limit. Since the above calculations are of the same order of magni- 
tude, it seems not improbable to conclude that the influence of thermal contrac- 
tion provided the conditions for the rotation of Mercury to be trapped into the 
3:2 resonance state. 
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